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Cosmic-ray (CR) protons can accumulate for cosmological times in clusters of galaxies. Their 
hadronic interactions with protons of the intra-cluster medium (ICM) generate secondary elec¬ 
trons, gamma-rays and high-energy neutrinos. In light of the high-energy neutrino events recently 
discovered by the IceCube observatory, we estimate the contribution from galaxy clusters to the 
diffuse gamma-ray and neutrino backgrounds. For the first time, we consistently take into ac¬ 
count the synchrotron emission generated by secondary electrons and require the clusters radio 
counts to be respected. For a choice of parameters respecting current constraints from radio to 
gamma-rays, and assuming a proton spectral index of —2, we find that hadronic interactions in 
clusters contribute by less than 10% to the IceCube flux, and much less to the total extragalactic 
gamma-ray background observed by Fermi. They account for less than 1% for spectral indexes 
< —2. The high-energy neutrino flux observed by IceCube can be reproduced without violating 
radio constraints only if a very hard (and speculative) spectral index > —2 is adopted. However, 
this scenario is in tension with the high-energy IceCube data, which seem to suggest a spectral 
energy distribution of the neutrino flux that decreases with the particle energy. We stress that our 
results are valid for all kind of sources injecting CR protons into the ICM, and that, while IceCube 
can test the most optimistic scenarios for spectral indexes > —2.2 by stacking few nearby massive 
objects, clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray 
and neutrino backgrounds in any realistic scenario. 
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1. Introduction 


The IceCube neutrino observatory at the South Pole has recently reported evidence of ex¬ 
traterrestrial neutrinos [jl]]. These neutrinos are compatible with a flux isotropically distributed in 
the sky, with astrophysical origin and with a possible cutoff at a few PeV. The origin of these events 
is unknown, however, the isotropic distribution in the sky of the observed events suggests that they 
might come from various extragalactic ~ 100 PeV cosmic-ray (CR) accelerators (see, e.g., [||]). 

As discussed in Ref. |Q], a multi-messenger connection between the measured neutrino fluxes 
and their photon counterparts could be crucial for unveiling the origin of the high-energy neutrinos. 
We here assume that the IceCube high-energy neutrinos have an extragalactic diffuse origin and are 
produced in proton-proton collisions. Therefore, they should be accompanied by gamma rays at 
a flux comparable to that of neutrinos. The extragalactic gamma-ray background (EGB) is the 
measured radiation that remains after subtracting all known sources from the observed gamma-ray 
sky, and it is due to the sum of contributions from different unresolved sources (see, e.g., 0]). If 
the IceCube neutrinos are indeed produced in proton-proton interactions, the same sources should 
also contribute to the EGB. 

Clusters of galaxies are the latest and largest structures to form in the Universe, where energies 
of the same order of magnitude as the gravitational binding energy, lO^^-lO^^ erg, are dissipated 
through structure-formation shocks and turbulence [|^. Even if only a small part of this energy goes 
into particle acceleration, clusters should host significant non-thermal emission (see, e.g., 0). The 
contribution of clusters of galaxies to the EGB has been discussed by several authors ||], and it 
has been argued that CR hadronic interactions, with the protons of the intra-cluster medium (ICM), 
in galaxy clusters could be responsible for a relevant fraction of the IceCube neutrinos 
However, such hadronic interactions could have a dramatic impact on the radio frequencies since 
secondary electrons are also produced in proton-proton interactions and radiate synchrotron emis¬ 
sion when interacting with the magnetic fields in clusters of galaxies (see, e.g., iHl)- Therefore, 


the radio emission from secondary electrons needs to respect radio counts of galaxy clusters [12]. 


In this proceeding, we summarise the main findings published in our recent article [ ]13[ ] where 
we estimated the possible contribution to the extragalactic gamma-ray and neutrino backgrounds 
from galaxy clusters assuming that gamma rays and neutrinos mainly originate in proton-proton 
interactions, while for the first time taking the constraints from the radio regime into account. 


2. Methods 

The total gamma-ray intensity from all galaxy clusters in the Universe at a given energy 
{AN ! AA At AE) is given by 


ly — 
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where the cluster mass is defined with respect to a density that is A = 500 times the critical 
density of the Universe at redshift z. Here, Ly is the gamma-ray luminosity of clusters, 14 is the 
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comoving volume, Dl(z) the luminosity distance, and d^u(M 5 oo,z)/dV'cd ^500 is the cluster mass 
function. The lower limit of the mass integration has been chosen to be Mjoo.iim = Mq, 

to account for large galaxy groups. The redshift integration goes from z\ = 0.01, where the closest 
galaxy clusters are located, up to Z 2 = 2. At the same time, the total number of detectable galaxy 
clusters in the radio regime at / = 1.4 GHz, above the flux limit F^in, is obtained as 


M.4(> 



dhi{FiA,z) dl4 
dVcdFi,4 dz 


dzdfi.4, 


( 2 . 2 ) 


where Fi ,4 = Li, 4 (l + z)/d7iD\^{z)'^, with L 1.4 the radio luminosity at 1.4 GHz. We compare this 
number with the radio counts from the National Radio Astronomy Observatory Very Large Array 
sky survey (NVSS) from Ref. [14]. The flux Fmin is defined as in equation (9) of [ [1^ by adopting 
a noise-level multiplier = 1 , which is appropriate, while slightly optimistic, for the low redshifts 
of the NVSS survey (0.44 < z < 0.2), and a typical radio half-light radius of /? 5 oo /4 [Hi- 

The function d^n(Fi. 4 ,z)/dV'cdFi 4 is obtained numerically from d^u(M 5 oo,z)/dVcdM 5 oo by 
calculating L\ a,{M^qq) from Ly{M 50 Q) as explained in the following. We introduce a phenomeno¬ 
logical gamma-ray luminosity-mass relation: 


log 10 


L^(lOOMeV) 

s-iGeV-i 


= Pi-hR2logio 



(2.3) 


where we omit the possible redshift-dependence for simplicity. The radio luminosity can be ob¬ 
tained from the gamma-ray one by using the “standard” formulae for the hadronic proton-proton 
interactions for which we adopt Ref. & considering a power law in momentum, p “p, for the 
spectral distribution of CR protons in clusters. In order to do so, we also assume that the magnetic 
field B in galaxy clusters is independent of the radius in the radio-emitting region. 

We note that the parameters Pi, P 2 , Op, B, and the fraction of loud clusters are degenerate 
when one tries to find the maximum allowed hadronic-induced emission. The most recent esti¬ 


mates suggest that the percentage of radio-emitting (loud) clusters is about 20-30% [|18[]. The 
subdivision of the cluster population into radio-loud and radio-quiet clusters is also reflected in 
the corresponding gamma-ray and neutrino fluxes. Therefore, from now on we refer to the two 
populations as “loud” and “quiet”, while we will mainly consider the overly optimistic case where 
all the clusters are loud (100% loud). In the following, to reduce the number of free parameters, 
we fix P 2 = 5/3 ~ 1.67; i.e., we assume that the hadronic-induced luminosity scales as the cluster 
thermal energy Pth Af^/Pyir where Pyir is the virial radius. The parameter Pi is then 

found requiring to not overshoot the radio counts from the NVSS survey and/or current gamma-ray 
upper limits, where for the latter we take the Coma cluster as reference case ||^. Eventually, we 
also calculate the total neutrino flux from all galaxy clusters in the Universe, similarly to what done 
in Eq. (^), as prescribed in Ref. [[T^]. 

We investigate different spectral indexes Op = 1.5, 2, 2.2, and 2.4, and different magnetic 
field values B » Pcmb>' P = 1 pG, and 0.5 pG. The first choice of the magnetic field can be 
regarded as conservative considering that the volume-averaged magnetic field of Coma, the best- 
studied cluster for Earaday rotation measurements, is about 2 pG [20]; while the last should be 


*ScMB — 3.2 X (1 -fz)^ qG is the equivalent CMB magnetic field strength such that is the cosmic mi¬ 

crowave background energy density. 
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considered optimistic with respect to current estimates. To clarify the meaning of the terms conser¬ 
vative/optimistic, note that the higher the magnetic field, the less room there is for protons, because 
the radio counts have to be respected, hence the lower the gamma-ray and neutrino fluxes. 


In addition to the above phenomenological method based on the luminosity-mass relation, we 
also adopt a more sophisticated approach for modelling the CR, ICM, and magnetic field distri¬ 
butions in galaxy clusters following Refs. |13, ^]. In this semi-analytical approach, the galaxy 
cluster population, as well as the ICM and magnetic field distributions, are modelled such that they 
correctly reproduce on average current observations. For the CR spatial and spectral distribution, 
we adopt the hadronic model proposed in Ref. |p^], which extends the semi-analytical model of 
Ref. [^] to account for possible CR transport phenomena. 


3. Results 


While we point the reader to our article [ [1^ ] for a more detailed account, we summarise here 
our main findings. Figure |I| shows the results obtained with the luminosity-mass relation approach 
for B = 1 /rG, which we consider the most realistic choice, for both the comparison of our models 
to the radio counts (on the left) and for the computed gamma-ray (in black) and neutrino intensities 
(in red) as functions of the energy (on the right), for the chosen values of Op assuming 100% 
loud clusters. For comparison, we plot the Fermi-Lwgt Area Telescope (LAT) data [21], and the 
four-years IceCube 1 a error band [||]. 




Figure 1: Total gamma-ray and neutrino intensities (right) due to hadronic interactions in galaxy clusters 
and the corresponding radio counts due to synchrotron emission from secondary electrons (left) for 100% 
loud clusters and B = \ pG obtained with the luminosity-mass relation approach. The neutrino intensity is 
meant for all flavours. All the plotted intensities respect NVSS radio counts and the gamma-ray upper limits 
on individual clusters. For Op = 1.5 and 2.4, the radio counts respecting the gamma-ray and neutrino limits, 
respectively, are below the y-scale range adopted for the panel on the left. 

For Op > 2, both the gamma-ray and the neutrino diffuse backgrounds are well below the 
Fermi and the IceCube data in all cases. For Op = 2, while the gamma-ray flux is always lower than 
the Fermi measurements, the neutrino diffuse background could represent a significant fraction of 
the flux measured by IceCube (for the cases of B = 1 ptG and 0.5 ftG). However, the case of 100% 
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loud clusters is not realistic. Therefore, in Figure we show the same as in Figure for Op = 2 
only, together with the more realistic case of 30% loud clusters. In the latter, galaxy clusters could 
make up at most about 10% of the neutrino flux measured by Ice Cube for B = 1 /iG. 



Figure 2: Same as Figure together with the case of 30% loud clusters for Op = 2. The remaining percent¬ 
age of 70% quiet clusters has been assumed to have Ly(lOOMeV) one order of magnitude lower than for the 
loud clusters. The 100% loud case is shown with lighter colours (i.e., in grey and orange). 


In the extreme case of Op = 1.5, we could explain the IceCube data by averaging over the 
corresponding energies for all cases, while respecting all other constraints from radio to gamma 
rays. However, we note that such a hard spectral index contradicts the most recent IceCube results 
suggesting a much softer spectral index [|^. 


Let us underline some of the assumptions made in obtaining these estimations. We do not 
assume any CR spectral cut-off or steepening possibly due to the high-energy protons that are no 
longer confined to the cluster [^, |^]. While this is not relevant when comparing with the Fermi 
data, it might be relevant for the high-energy neutrino flux and, therefore, our results should be 
considered conservative in this sense. Additionally, we omit the absorption of high-energy gamma 
rays due to the interaction with the extragalactic background light because this becomes relevant 
at energies above few hundred GeV for low redshift sources such as galaxy clusters [24\. Another 
important assumption, worth discussing here, is the slope of the luminosity-mass relation P 2 = 5/3. 
While this choice is appropriate for CR protons injected by structure formation shocks, it could be 
different for other CR sources. We, however, estimate that a much steeper slope would imply lower 
gamma-ray and neutrino fluxes, as it would require to lower Pi to not overshoot radio counts, and 
that a much flatter slope could realistically boost the total extragalactic neutrino flux to be at most 
about 30% of the IceCube flux for Op = 2. However, a luminosity-mass function with a slope both 
much steeper and much flatter than 5/3 would contradict current knowledge of the scaling relations 
of diffuse radio emission in galaxy clusters 1 1^]. 


We can conclude this section saying that amongst all the considered cases that respect both 
radio counts and current gamma-ray upper limits, hadronic interactions in galaxy clusters can 
realistically contribute at most up to 10% of the total extragalactic neutrino background, while 
contributing less than a few percentage points to the total extragalactic gamma-ray background. 
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Moreover, we do know that not all the observed radio emission in clusters has a hadronic ori¬ 
gin |]^, ^], therefore, the simplified requirement of not overshooting the NVSS radio counts on 
clusters leads to optimistic results. This implies that even our results, which respect both NVSS 
counts and gamma-ray limits, should still be considered rather optimistic. We note that, owing 
to our simplified approach using a gamma-ray luminosify-mass relafion, fhese conclusions can be 
generalised fo any source of CR profons where fhese mix and hadronically interacf wifh fhe ICM 
of galaxy clusfers, such as fhose injecfed by strucfure formation shocks and AGNs. In facl, for 
any considered source of profons, fhe resulfing secondary emission musf respecf bofh radio and 
gamma-ray constraints. 


4. Future detection prospects 


Recently, Ref. [23] have provided neutrino flux upper limifs on a sfacked sample of nearby 
galaxy clusfers (Virgo, Cenfaurus, Perseus, Coma, and Ophiuchus) following predictions provided 
by Ref. [^]. These upper limifs are /v,ul( 250 TeV) = 6.9 x 10^^° and 7.7 x 10^^° cm^^ s^^ GeV^' 
for fhe cases where CR profons are supposed fo be uniformly disfribufed wifhin fhe clusfer virial 
radius (A) and where CRs are assumed fo be disfribufed like fhe ICM in clusfers (B), respecfively. 

The above neufrino flux upper limifs were obfained by assuming a specfral index of « 2.15, 
so we can compare wifh our resulfs for Op = 2.2. For fhe slacked sample of Virgo, Cenfaurus, 
Perseus, Coma, and Ophiuchus, we eslimale lhal fhe upper limifs from Ref. are jusf a factor of 
1.3 (A) fo 1.5 (B) above fhe maximum allowed hadronic-induced neufrino flux. When Op = 2.4, fhe 
maximum allowed flux for fhe slacked sample is one order of magnilude lower, while for Op = 2 
if is one order of magnilude higher, wifh respecf fo Op = 2.2. While we slress lhaf any realistic 
modelling of fhese objecls should consider fheir specific characferislic, in particular fhe presence 
or lack of diffuse radio emission, we can conclude lhaf IceCube should be able fo pul conslrainfs 
on our mosf oplimisfic scenario wifh Op = 2 and on fhe Op = 2.2 case in fhe very near fulure, while 
fhe case wifh Op = 2.4 will much harder fo fesf. 


5. Anisotropies in the gamma-ray background 


The aulhors of Ref. pSQ have shown lhal a major fraction of fhe EGB is made by discrete 
sources and fhe measured level of anisolropies is consislenl wifh prediclions for gamma-ray blazars. 
While we showed lhal clusfers are nol fhe dominanl conlribulors to fhe fhe lofal gamma-ray and 
neufrino fluxes, Ihey could make subslanlial conlribulions to fhe EGB anisolropies. 

Eigure ^ shows fhe angular power speclrum eslimaled using our semi-analylical model and 
compared wifh fhe EGB measuremenf, and upper limifs obfained once fhe main blazar componenl 


is sublracled []28[]. Our prediclion is abouf one order of magnilude less lhan fhe Eemri-EAT upper 
limil. We underline lhal in scenarios where fhe lofal galaxy clusfer intensify is much higher lhan 
for our semi-analylical model, as is potentially realised for some of our simple phenomenological 
models based on fhe luminosify-mass relation, fhe angular power speclrum could be a powerful 
discriminalor, as powerful as radio counfs. 
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Figure 3: Gamma-ray angular power spectrum for emission resulting from proton-proton interactions in 
galaxy clusters in the energy range 1 — 10 GeV. We show the result for our semi-analytical approach for 20% 
and 40% loud clusters. We also show both the measured EGB anisotropy and the upper limits obtained once 
the blazar component is subtracted [^. 


6. Conclusions 

We estimated the contribution from hadronic proton-proton interactions in galaxy clusters to 
the total extragalactic gamma-ray and neutrino fluxes, while including radio constraints for the first 
time. We modelled the cluster population by means of their mass function, and make use of a 
phenomenological luminosity-mass relation applied to all clusters constructed by requiring both 
radio counts and current gamma-ray upper limits to be respected. We also adopted a more refined 
approach that employs a semi-analytical model for CR, ICM and magnetic field distributions in 
clusters. With these we showed that galaxy clusters can contribute at most up to 10% to the neutrino 
flux measured by IceCube, while contributing much less to the EGB. Nevertheless, IceCube can test 
the most optimistic scenarios for spectral indexes > —2.2 by stacking few nearby massive objects. 
We also discussed that galaxy clusters could substantially contribute to the EGB anisotropy because 
they are fewer in number than other astrophysical sources and, therefore, are expected to be more 
anisotropic. Concluding, our results put earlier work, which turned out to be overly optimistic in 
estimating the galaxy cluster contribution, into prospective. 
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